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ABSffiACT 


In this V7ork, thermal fatigue analysis of pressure 
vessels with hemispherical, elliptical, torispherical 
and toriconical heads has been studied using finite ele- 
ment method. Inside of the vessel is subjected to pre- 
ssure and cylic temperature variation of heating, hold- 
ing and cooling; and outside is insulated. V/ei^t of 
the vessel has been considered and isoparametric eight 
noded quadrilateral ring elemcsits have been used. 

Stress variation is studied in all the vessels. !.as 
function of heating/cooling rate and hold time. It is 
observed that maximum stress variation occurs in the 
skirt region near the head and is tie dominant factor 
for the thermal fatigue life of the vessels with all 
heads. Vessel with hemispherical head is the best design, 
and others in this order are elliptical, torispherical and 
toriconical heads. 

finally the design cmnges are considered to reduce 
the stress variation in the skirt. One is to change the 
position of the skirt and other is to provide a portion 
of it a temperature variation as that of the inside of 
the vessel. The second method reduces the stress variation 
considerably. 



CIL'LPiTS’ I 
III2R0DUGTI0N 

In recent years of research and develoiraent considerable 
attention has been focused on tJie design of pressure vessels 
because these are the means by v/hich energy systems controlled, 
chemical and petroleum processes operated. Modern chemical 
processirjg plants and power generation industry demand the use 
of very high pressures and temperatuces. Many chemical pro- 
cessing press-ure vessels are subjected to a large number of 

I 

temperature cycles, hence to stress c 3 )cleso Ultimately the | 

i 

vessel is liable to fail 'by thermal fatigue. Thus in such 

1 

cases tbe designer is charged v/ith the task of doing thermal 
fatigue analysis. This can only be done by a careful stress j 
analysis of the entire vessel. 

1,1 HIMIOUS VJORK 

Pressure vessels and its components are designed based 
on the standard codes like ASJ-iE, i3S §tc. However, the design- 
er is often faced with complex loadings which are not covered 
by the cades. Thus certain designs fall beyond the scope of 
the code. One such loading is that of the inside surface of 
the vessel being subjected to cyclic temperature variation. 

This will induce transient thermal stresses vdiich depend on 
the temperature attained and -Hie thermal gradients caused. It 
is not possible to estimate the temperatures and the corres- 
ponding thermal, stresses, as closed form relations are not 


1 
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available. Go a design of tbe vessel is completed on the basis 
of mechanical loads. 

When once this design is completed, a detailed thermal 
stress analysis is needed when the vessel is subjepted to 
thermal transients. In thermal fatigue problems, the designer 
is interested in loiov/ing what is the critical section, where the 
stresses are going to bo ma:cimum. Also for calculating the 
allowable number of thermal cycles, knowledge of maximum alter- 
'nating stress or strain amplitude is essential. 

S'bress analysis of pressure vessels can be performed by 
the follov/ing methods: 

1 « Analytical Kethods 

2. Experimental Kethods v/hich include 

a. Strain gauge measurement 

b. Photoelastic method 

c. brittle coating method 

3. numerical Techniques. 

Analytical methods require a general solution of the equi- 
librium equations. But sometimes the problem becomes too complex' 
and beyond analytical solution. Then one has to resolve to experi- 
menjbal methods. But the ojqjerimental meidiods pose a lot of diffi- 
culties in the fabrication of the model and creating the required 

of 

type of loading. This needs a lot/skill and is highly costly in 
time and money. 

Thus in recent times, numerical Techniques such as finite 
element melhod has been used for the stress analysis* In E®!, 
the continuum is divided into a Unite number of elenents. A 
solution for the displacements is assumed initially in the form 
of a polynomial without giving any consideration to the boundary 
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conditions of Ihe problem, ITov; tising eithgL- Tariational princi- 
ples or v/eighted residual methods, typical finite element equa- 
tions between stiffness, displacements and load mati‘ices are 
formed, These elemental equations ai’e suitably assembled and 
global equations are formed. Then by applying the boundary 
conditions, the system of equations .are solved for the displace- 
ments and Icter strains and stresses are found at various points. 
In this way a general 5'2I-I program can handle a -mrie-ty of pro- 
blems of Solid Hechanics, Heat Transfer and PI uid Mechanics , 

The development of modern computers v^ith tremendous speeds of 
calculation made PEM to emerge as a powerful numerical technique 
in the field of stress analysis. The results from the PB-I are in 
excellent agreement v;ith theoretical exact solutions available 
for some pr obi cans, 

hanger [|l]] reported thit since most pressure vessels are 
subjected to limited number of pressure and temperature cycles 
during ihc life time, considerable design effort could be saved 
by defining the conditions which do or do not requrre a fatigue 
evaluation. The author in his studj'' concluded that a fatigue 
evaluation is needed if during normal operation, 

1, Pressure fluctuates tlirough a range exceeding 
t\^^enty-five percent of the desigjn pressure, ca: 

2. The temperature difference between any two adjacent 
points in the vessel changes by more than 22®G (40®P), 

This malces it possible to eliminate a great deal of detailed 
calculations for vessels, that ore subjected to pressure and tem- 
perature cycles of moderate severity. 
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Pickett and Grigory [2] reviewed the fatigue evaluation 
procedures and the fracture mechanics approach to fatigue life 

analysis, They made important remarks on crack initiation, 
crack propagation £ind about the fatigue test specimens. 

They mentioned tliat specimens loaded in bending are better 
compared with axial or push pull type specimens. They con- 
cluded that the surface strains at the test section on bending 
specimens shake down to stable values after a, fev/ cycles and 
, can be measured easily. 

Host studies of lov/ cycle fatigue behaviour have been 
concerned with evaluating the effects of cyclic frequency and 
holding time at peak strain or stress, 

Conv/ay, Berling and Stens £^2 an experimental work 
using servo controlled, hydraulically actuated fatigue machines 
v/hicdi subjected inductively heated specimens to a push pull 
type of loading. They developed a diametral strain extenso- 
meter and the signal from ihis device and from the load cell 
in series with the specimen was' fed to a specially designed 
computer network to provide for axial strain control. 

They reported that hold period at peak stress or peak 
strain leads to reduction in the cyclic fatigue life. They 
also reported that 1he time to fracture is proportional to the 
length of the hold period. This effect is just opposite to the 
number of cycles to . • fraotui'e. As the length of the hold 

period increases, number of cycles to cause failure decreases. 

Saxaph, Kushwaha and ICakodkar [4] worked on transient 
thermal stress analysis of pressure vessels using finite element- 
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method o They studied the problem of thermal shock, one of 
the most serere loadings experienced by light water cooled 
reactor vessel, when it is subjected to emergency coolingo 
They modelled the vessel as an axisymmetric geomelacy with 
eight noded isoparametric finite elements. The transient 
temperatuces were calculated by Crank-ITicholson* s scheme. 
Detailed temperature and stress plots were given. The 
authors had taken a cylindrical vessel with a hemispherical 
head which is the moot ' one in nuclear industry, 

1.2 PRESSm^ v\DRK 

In design of vessels subjected to mechanical (presstire) 
and thermal loads, both should be {^ven proper importance. 
Because, the principles upon which -aieimal stresses can be 
mitigated are many times different to those involved in 
designing for mechanical load, Bor example, making of 
sections of large thickness can reduce the operating stresses 
for mechanical loads. But this approach may increase thermal 
stresses, as hi^ tliormal gradients are set up in relatively 
thick vessel walls. 

In chemical and petrochemical industries, pressiare 
vessels with hemispherical , elliptical, torispherical and 
tor iconic al heads are q.uite common,' Hence in the present 
work, thermal fatigue analysis of piressure vessels with four 
types of heads mentioned above has been studied using the BBM. 
Thermal stress calculations need prior knowledge of tempera- 
ture distribution. Both temperatures and thermal stresses 
have been determined by the finite element method, Mlhougii 
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still it is more popiiLar to iise finite difference method for 
temperature calculations hut here I’M has been usedo As ¥M. 
is ■very straight forward method and implementation of boun- 
dary conditions is an extremely simple affair o, In addition 
use of PM for temperature and thermal stress calculations 
reduces •the overall data preparation as same mesh is used 
for both* 

It is known that hi£^er order’ elements give •the same 
accuracy with fev/er number of clQiicntc, This also requires 
much less data preparation. In addition pressure vessels 
have many curved surfaces. So the isoparametric eight noded - 
quadrilateral ring elements have been used for temperature 
and thermal stress calculations. 

In the present work, preccvire vessels are supported 
by skirts. Insides of the veci’els are subjected to cyclic 
temperature variation of heating, ‘ holding and cooling and 
outsides are insulated. 

Allowable number of cycles as function of heating/ 
cooling rate and hold time are calculated. It is found 
that maximum stress variation occurs in the skirt near the 
head and is the controlling factor for the life of the vessels. 
So design modifications are considered to reduce these stress 
variations in the skirt and thx-is improve the life of the 
vessel. 



•prs .j 

'Vi/ ^ i Jk ■I•JL ^ 1 ^ aia, 

l-GlL.UL^/JICil 


Finite element method ic tieed for tlie tlierrnal f?>,tiguo 
analysis of prercnrc v"GSSv.;le ir. tliic \rorko GalciiLa.tion of 
the thermal stresses reiiiiircs the tempera time distrihutiono 
Fii'st the method to dete-vine the tempera txiro disti’ibution 
is civen and then the method to determine thermal c tresses 
is given. Details of these methods ai*e avo-ilablc in the 
te 2 :t [ojbut are iiiclndod he..-c to make the thesis self- 
contained. 

As the ain of ihu thesis J to study the tlibrnal fatigue 
life of preocuro vessels v/iLh various heads, desi^in of vessels 
with elliptical, torispherical and toriconical heads is also 
given in this- chapter* 


2.1 IRAITSIiadT gZ'.l.fUfURJ C;A.3ULl!fI0h 

Pressure vessels under consideration are axisymmetrio 
bodies subjected to ao:isinmiietric temperature distribution on 
the inside and incvilated on the outside* (Thus the temperature 
distribution is indcpcndexit of © coordinate* Therefore, the 
governang differential equation for the temperature T is [j[] 

1 1 




Sr 


■* r ^r 


a «>t 


( 1 ) 


v/here a is the’ thonii'il dijdhuaLvlti'. 

For this problem f ini Lc c=>] cinart.'s become ring elomonts. 
To take cai'e of cin.‘ved bOkUid -n-ies, isoparametric eiglxt noded 
quadrilateral , el orients arc tidfcii as cliovni in Fig* 1, 
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[temperature distribution over the element is taicen as 

= LIJ (2) 

where are the interpolating' functions and are 

the temperatures at the nodes. 

Substitution of 3qn, (2) in Eq.no (1) £;ives the residue 
of the element, fhis residue is minimised by the G-alerkin 
method and finally one ceis 

[0]^®> + [K]^®^ fljCne) = (3) 

■ where 

= /J's * WL:tJ to da 

= //rC{i;,rJ + jH.zj (1!,^) dr dz 

* 

and 

= /.{i:j(^)rdh ■ (4) 

[c] is the thermal capacitance matoix, is the 

thermal stiffness matrix and is the boundary heat 

flux matrix.- 

2.2 fllERIhH S!I!aB5SES. mi i OP J lULn^ICh 
The potential to bo minimised is 

1 

I = i-J«r6 dv - - /Sj^U^ds C5) 

Here the first term represents the strain ener^jy* the second 
term represents the v/orh done by bodj forces and the "feird 
term represents the'v/orl: dcaio surface foi'ces, 

fhe pressure vessels under consideration are axisymmetric 
bodies subjected to axis^TKiiCtric pressure and temperature loadings, 
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Thus the displacenents, Btrains and stresses axe independent 
of Q coordinate', i Thus for tjiis problen also finite elements 
become ring elements as for temperature distribution problem. 
Same isopareoaetric eight noded quadrilateral elements shoivn 
in Pigo 1 are taken for this problem, but here each node will 
have the degrees of freedom u and v. 

The displacements over the fijiite element are taken as 


0 


c o o o 


0 •^'*-1 ^ 0*0 ooo 

I 


£uj ^ ® ^ = [U (r , z )] ^ 


Here u is the displacement in t3ic r-direction and v is tho 

j 

■displacement in the s-direction, 

The strain matrim for this, axic^^maetric problem irJtC63 


W= 


6" ' - 


Sz ^r 


Using Eqn, (6), the strain matrix {€} , Eqn* (7) becomes 


= [D] = [jiji 1;^}, 
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v/here 


^"i 

0 

1 i.; 

r ‘ 


0 

0 s 

OH. 

1 


Os 

5r" 


( 9 ) 


The stress-strain ms. trices are related to ori,Q:i other as, [^63, 


vdioro 



is the olaoticity niotri:: and. is the 


(10) 

thermal * 


strain matri::, and these are , 

Ti 



S(1 - « ) 

TT+ o)(r-Tn 


■0 ^ 

t :::? 


1 


T35 


0 

0 


sjmi 1 0 


1-20 


oc T 




T 

0^ T 


( 11 ) 


( 12 ) 


L °J 

Here is the coefficient of linear thermsil expansion and T 
is the temperature rise * 
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Using Eq,ns, (6), (8) and (10), Eqn, (5) becomes, 

I = ffiBf [E] [Blrdrdz 

- 27r^[33“ [D] {€Q}rdrdz 

“ prds 

- (iij^^®^ I’rdrdz (13) 

By applying the Rayleigh. Ritz method, one gets. 


where 

[e]'®? =2nff'CB]’^ M [3] rdrdz 

{Eq}^''®^= [d]{€q} rdrdz 

27T^Cn 1^ prdB 

2^11] ^ Erdrdz (15) 

Here [^E]^®^ is the stiffness matrix, is iiie thermal 

load matrix, consistent surface load matrix and 

is bodj’' force mati’ix. 

In this present problem surface loads are the internal 
pressures acting normal to the surface as shown in Pig, 19o 
The consistent surface load matrix is , 

{p}(^®) =27TfCHj^ {pfj’''^® (16) 

If the pressure is acting along the edge *5^ » + 1 , then ' 

% ■“ P dB ~ ^ dB 

3s 

or p^ dB s= pdz *= P “ "* = *• P 

Substituting Eq.n. (17)-inEqno (16), one gets. 



( 17 ) 
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{p}(ne) 





r d 





(18) 


Similar expressions for other edges' is = -i and ^ = + 1 
can he writteno Value of this integral is determined hy Gauss 
Legendre numerical integi'ationo 

Body force in the present problem is gravitational 
force. Thus the body force matrix becomes, 



(19) 




2.3 


^ /XL ^ 

^ V> M ■ S.I I 


2,.3o1 Determination of Matrices 

All the matrices developed in Sections 2,1 -and 2.2 are 
evaluated by Gauss Legendre numerical integration scheme. In 
this scheme the quadrilateral is first transformed into a 
sqtiare of side 2. Then the ftuiction inside the integral is 
evaluated at each Gauss point, and multiplied by the corres- 
ponding weight. Sum of these products give the value of the 
integrals and thus the matrices. 2x2 point integration has 
been used in this work. 

Determination of thermal load matrix {I'q} of Eqn. 

(15) needs some comments. In this matrix, one needs the 
temperature rise T. Reference [6] suggests that this be 
taken as i-average temperature rise of the nodes of the element. 
This does not seem to be an accurate method for evaluation of 
this matrix. It is more accurate to determine the temperature 
rise at every Gauss pointf; and use it while doing the numeri'tjal 
integration. 
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2o3o2 TeimDerature Determination 

Eq.no (3) are first order simiiltaneous differential 
equations. They are solved here by finite difference 
method. Thus T^ at any time at i node can be 'expressed 
as , 

mP+1 _ mP 

m _ _i i 

" At 

and 

T^ = 0 + (1-0)T^ (20) 

' . 

where p refers to previous time step and p + 1 refers to 
the current time. 

Using Eqn. (20), Eqns (3) become. 


(-4 + 0 [K]^®^) {t} 


(e) 


,(ne) 


P+1 = (^ - (i.0)Iic]^®)){t} 

+ 


or 




(ne) jp . 

(21) 

l(ne)|p 

(22) 


These matrices m and [iTjcan be evaluated once © and At 
are fixed. Those equations are assembled for the whole domain 
and boundary conditions are applied. The assembled equations 
are simultaneous algebraic equation and are solved by Gauss 
elimination method. This way temperatures at all the time 
steps of the cycle are determined. 

Value, of 0 was talcen as 1,0 in this work and is discussed 


(ne) 


in Chapter III, 



14 


2.3*3 thermal Stress Determination 

Pinitc element eqnco (14) are aosembled in the usiial 
way and hoimdary conditions applied* (These are simnlta~ 
neous aloehraic equations and are solved by G-aussian eli-. 
mination for displacements. Stresses are determined using 
Eqn. (10). 


2.4 DESIGN 0? VARIOUS PilSSSURE YBSSBL HEADS 

A design of the pressure vessel with hemispherical 

heads done by Bharat Heavy Plate & Vessels has been available 

with Department of Hechanical Engineering, Indian Institute 

of (Technology, kanpirc. Elliptical, toriopherical and tori- 

conical heads are designed for the vessel using ASME code 

given b;y Chuuo & Eber [11 j|. 

Design details available for vessel with 
hemispherical heads 


Design pressure 

"forking pressure 

Code followed 

Material of the vessel 

Design stress intensity of 
the matoi’ial 


28 MPa 
26 lyiPa 

ASME 1980-81, SEC VIII 
SA 336 Gr E11 
155 MPa 


9 o 4936x10“^ mVsec. 


12.474x10' 
7860 kg/m' 


-6 


Yield strength 275 MPa 

Ultimate tensile strength 480 MPa 

(Thermal diffusivity 
Coefficient of linear expansion 
Density of the material 
(The Pig. 2 gives the other design details of this vessel 
with hemispherical head,. The thickness of the hemispherical 
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head has been calculated by the following relation {il^ 


t = 


PR 

2SE -'0.2? 


(25) 


where P = internal prccs'ure 

R = inside ra.dius of head as shown in Pig. 2 ’ 

S = allowable stress 
E = joint efficiency. 

Elliptical head: Thiclmess of the elliptical head is 
• given by [11] 


t = 


PD 

2SE - 0o2P 


(24) 


v/here D is the head skirt diameter. This head has been 
designed for a mjor axis to minor axis ratio of 2:1. Thick- 
ness comes out to be I^Omm. Other details are shovm in Pig, 3 » 
Torisphorical head : The thickness of the torispherical ■ 
head is given by [11] 


0.885 PI 
“ SE - 0.1P 


(25) . 


where I is the inside crovm radius, I has been taken 974mm 
in this design. V/ith tliis, head thickness comes out to 155mm. 
Knuckle radius has been taken as 6^ of the inside crown radius 
as re(iuired by ASME code [ll] . Other details are given in 
Pig. 4. 


Toriconlcfil head : The thickness of the toriconioal head 
is given by [ll] 

PD. 

^ ( 26 ) 


2cos «< (SE-0,6P) 

where D-j is the inside ’cone diameter at the point 'of jfeangency 
to the laiucklei is one half of the apex angle. In this 
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design has been talcen as I455iiiiii and «C as 50®. Inside 
knncklc radius has been talcen as 65^ of the outside diameter 
of head skirt* Other details are given in Pig. 5 o 



CIL4PTSII III 


RESULTS and LlSCUSSICir 

As already mentioned, finite element method has heen 
used for the study of thermal fatigue life of pressure 
vessels \vith different types of heads. The stress states 
in various heads have been studied. The effect of rate 
of change of temperature and hold period on the thermal 
fatigue life of the pressure vessels has also been studied, 

5.1 OHECEIITG THE VALIDITY OF TIGS HlOGR-al^ DEVELOPEL 

To check the validity of the finite element program 
developed, a test problem given in £53 was used. Also the 
problem of thick sphere under internal J)ressure and steady 
state temperature given in \'ja.s used* A good agreement 
of the results • with the a.vail^le solutions was observed, 

5.2 dRAESMT TH^IPERATURE DISTRIBUTION 

Eor transient temperature calculations, both explicit 
(© s 0) and implicit (Q 0) were tried v/ith, to find the 
value of Q which will give minimum or no oscillations of 
temperatures. Generally it is accepted that © greater than 
1/2 t^ill give less oscillations £63* Values of © as 0, 
(Galerkin’s scheme) and 1 (fully implicit) were tried. 

It was found that fully implicit scheme of © = 1 gave 
almost zero oscillations. So in this v/ork this scheme has 
been used through out. 
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SOKESS C O MPUTAlEIOir i^D ASSESSMENO? 0? PAgIGUE LIPE 

After calculating tbe temperatures in the vessel at 
various time steps of the thermal cycle, stresses v/ere 
determinedo The temperature and stresses in the vessel 
were determined using 2x2 point integration* It was 
observed that the stresses in the vessel reached an extreme 
at the end of the heating portion of the cycle and reached 
another extreme at the end of the cooling portion of the 
cycle* This was found to be true for all types of thermal 
cycles i*e,, with or \^dthout any hold period* As only 
maximum alterating stress amplitude v/as needed for the 
fatigue life evaluation, stresses were calculated only at 
these two points of the thermal cycle i*e,, at the end of 
heating and cooling portions* In this way, a lot of computer 
time was saved by avoiding the stress calculations at every 
instant of the thermal c^’cle. 

Further it was observed that the principal stress . 
directions at the above extreme points of thermal cycle were 
different. So to find the stress amplitude, the method 
given by ASME code [12] vra.c used, which is as follow: 

The stresses trz extreme 

points are calculated and subtracted* These stress diffe- 

" I 

rences are denoted by Ve** <rz* prin- 

cipal stresses corresponding to these stress differences 
are denoted by «r2» tf-j their differences as 

<ri* - ^"2* “ ^3* “ 

maxamtm difference gives the maximum stress ' amplitude* The 
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allov/alDle number of fatictie cyclec are dcberiained using the 
design 'fatigue c\u‘ve given in ASIiE code £123* 

5o4 SIERESS S^aTE 111 YiRIOL'S VESSELS 

All vessels are subjected to an internal pressure of 
26 MPa* Insides of the vessels are subjected to tempera- 
ture cycles, typical one is showi in the Fig, 6, and out- 
sides are insulated* Temperatures vary from 5®C to 245®0 
on the inside* Two rates of heating/cooling of 1®C/m, and 
1.5®C/m are studied. For each heating/cooling rate five 
hold periods of 0,1, 2, 3 and 4 hours are considered. The 
thermal cycle v/ith licating/cooling rate of 1.0®C/min and 2 
houi's of hold period is aalled here as a typical cycle, 
3*4.1 Stress Distributiorj hue to Mechanical loads 

Maximum' principal stress in upper head, cylindrical 
portion and in lov/er head is shov/n in Table 2. 

Table 2: Maximum Principal Stress, MPa 


Vessel v/ith 

Upper Head 

Lower Head 

Cylindrical 

Hemispherical 

Head 

152 

173 

136 

Elliptical 

Head 

184 

182 

142 

Torispherioal 

Head 

175 

189 

141 

Toriconical 

Head 

117 

126' 

141 
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It is seen that the principal stresses are norc in the 
lov;er head than in tlie upper head except in the vessel v/ith 
elliptical heads o These stresses are essentially same in 
the cylindrical portion. Principal stx’ess in upper head 
is Eiaximtim in the elliptical head and niiniinum in toriconi- 
cal head. Principal stress in the Mov/sr head is maximum 
in the torispherical head and minimum in toriconica.! head. 

It is also observed that this maximum principal stress 
occured on the inner side for all the heads. 

Maximum principal stress across AA in the upper heads 
is shov.Ti in Pig. 7. It is seen •tho,t it essentially remains 
constant across the thichnass in hemispherical head, and 
fall very sharply in the other heads. 

5 . 4.2 Tonperature Distribution 

Tempera tuire distribution in the vessel with spherical 
heads ‘at the oiid of hcati.ij is shovn in Pig. 8 and at the end 
of cooling in Pig. 9. These are for a typical thermal cycle 
in the similar ’./ay tempera t’ures are calculated lor other 
vessels and subsequently used for stress calculations. 

5.4*3 Stress Distribution Due to Fiechanical and Thermal Loading 
Maximum principal stresses in upper head, cylindrical 
portion, lov/er head and skirt at the end of heating and 

cooling arc given in Table 5. These are for typical thermal 

✓ 

cycle • 

It is observed that generally maximum principal stress 
occured in lov/er heads both at the end of heating and cooling. 
Also it is obsorved that iK-riiiruUU principal stress oocurBd at 
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the end of cooling in all .oections of the vessles except in 
the skirt. It is seen that the naximmi princixjoJ. stress in 
the cylindrical ijortionc of al.1 the vessels at the end of 
hoatina roiiiainod noro ir c..-jac . Thv., raaic is- ohser-vod 

at xh: end of co-.li]:^^, Jx i.< ^ be .noted thrxt maximnin 
values at ifche end of heatin/i and cooling many a times do not 
occur at the same point, 

Maximim principal stress across AA in the ujaper heads 
is shovm in Pigs, 10 and 11. It is obooi’ved tliat maximum 
principal stress in all the heads occxirjed at the inside 
surface, and fal'tc ch.irply across the thickness. 

Maximum principal stress across the cylindrical portion 
BB of the vessel v.lth sphoric.'-l head is shovm in Pig, 12, 

In Pigs. 13 and 14 ohe variation of maximum principal 
stress in varior.s heads ii a topical thermal cycle is shown, 
These can be referred to as tjpicol stress cycles produced 
in various heads and are draxn at a point' in an olaaent 
v;hcre stress arajl.i tr<' c j;: cbrerx’'ed to be maximum. It is 
observed fro:: these ?ijr_.*:r tliat the hemisphErical head is 
having the least stress amplitude. It is observed that the 
stress amplitude is increasing in the order of elliptical, 
torispherical and toriconical heads. In Pig, 15 the varia- ■ 
tion of maximum principal, stress in a typical thermal cycle 
in the cylindrical portion of the vessel v/ith hemispherical 
heads is shovm. 


♦ i 

a.**! 


UlL 

Ip - 


aPjm 
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It is observed from thse graphs that dtiring the heating 
portion of the cycle, "ihe thermal stresses in the element 
considered are compressive and dircing the cooling portion 
they change to tensile stresses# 

5o5 THERI'IAL FATIGT3S LIFE 

In this study, two heating/cooling rates of 1.0 and 
1.5®G/m v/ith five different hold periods as detailed in 
Table 1 are considered. Thermal stresses and maximum stress 
amplitudes are calculated as explained in Section 3.5o 

Maximum stress amplitudes for these thermal cycles 
for vessels v/ith different heads are given in Tables 4 
and 5 for 1 and 1.5®C/m respectivoLy. It is seen that these 
amplitudes hardly change in the upper head and cylindrical 
portions, slightly change in tlie lover heads and considerably 
change in the skirt regions, for all the vessels. Also it 
is seen that these amplitudes are more in the lower heads 
than those in ihe upper heads due to the proximity of skirt 
near the lower heads. It is observed that maximum stress 
amplitude always occurs in the skirt region near the vessel. 
Values of these amplitudes are shown in Fig. 16 for Tarious 
vessels for different hold periods. As expected the maximum 
stress amplitudes increase with increase in hold period and 
heating/cooling rate. It is seen that these amplitudes in- 
crease in the order of hemispherical, elliptical, tori- 
spherical and toriconica.1 head vessels. 
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Fati^jtie life of various vessels i»e«, allov/able number 
of fatigue cycles is calculated using tbe design fatigue curve 
given'in ASME code 0^1*- Fatigue lives of various vessels for 
different beating/cooling rates and hold periods are shown in 
Table 6« For heating/cooling rate of 1®C/ni these values are 
also given in Pig« 17 as functions of total cycle timeo For 
heating/cooling rate of 1 and 1,5°C/m these values are given 
in Figo 18 as functions of hold period. For vessels with 
hemispherical heads, with heating/cooling rate of 1*0®C/min, 
the fatigue life decreased by 23'^ from a no hold period to a 
hold period of fouc hours. For the same conditions mentioned 
above, the vessels with el3.iptical, torispherical and' tori- 
conical heads ohov/ed a decrease of 28^, 345^ and 365 ^ respectively. 

With a heating/cooling rate of 1,5®C/inin, a change in 
hold period from zero to four hours, vessels with hemi- 
spherical, elliptical, torispherical and toriconical heads 
show a decrease in fatigue 1 ife of 23^, 29%, 3A% and A^% 
respectively. 

For a hold period of 2 hours, cliange of heating/cooling 
rate from 1.0 to 1o5®^j/m vessels with hemispherical, elliptical, 
torispherical and toriconical heads shovi a decrease in the 
fatigue life by 50?^, 505^, 26%, 14^ respectively. 

3.6 ICBTHOFG TO FvBDUOE THE STRESS AT^llPLITUPE 11 TliB SKIRT REGION 

Throughout the analysis, it is observed that high stress 
amjjlitudcs are developed in the skirt region adiacent to the 


2'4 


Tiea-do I'lie cxictin^ decii^'n oj? tho r.kirt caused large thermal 
gradients at the alcirt vessel junction, thereby producing 
high stress amplitudes, ibus in the entire analysis the 
skirt is found to be the predominant factor in reducing the 
fatigue life of the vessel, because of the skirt, higher 
stress amplitudes arc produced in the bottom head also. So 
efforts are made to reduce the stress a;.iplitude in the skirt, 
thereby increasing the fatigue life of the vessel, Ihe tv/o 
methods which successfully reduced the stress amplitude in 
the skirt are as belov/t 

1 , Changing the position of the skirt, 

2. Providing a suitable temperature distribution 
in the skirt, , 

In the existing designs the skirt is attached at the 
head-cylinder junctions. In the nev; design, the position of 
the skirt is changed caid ' attached only to the head region 
as shown in the Pig, 2. This change in design is studied 
in the vessel with hemispherical head only. Per the no holid 
time thermal cycle with 1,0®G/min heating/cooling rate, the 
nev/ design reduces the stress amplitude by about lOjS, But 
even this 10^ decrease in stress amplitude does not increase 
the fatigue life much. So a second method mentioned above 
has been tided with. 

In this second approach, a favourable temperature dis- 
tribution is provided at the skirt to decrease the existing 
high thermal gradients'. Hence one fourth length of skirt 
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adjacent to the head has Deon subjected to the same thermal 
cycle as that of the iJisido of the pressure vessel. The 
required tcmpex'atni'e variation v/ith the desired rate of 
chan:;:e of temperature maj’’ be brought by using electric 
heating coils around that skirt. This change considerably 
decreases the thermal gradients at the skirt vessel junction 
and hence the stress amplitude. For a heating/cooling rate 
of 1®C/m, maximum stress amplitudes in upper head, cylinder, 
lower head and skirt of various vessels are shovm in Table 7. 
Results for these cases v/ithoxit skirt heating are also given 
in this Table. Maximum sti-ess amplitudes decrease considerably 
in lower heads and skirts. Decrease in skirts is almost 
dramatic. For vessels vdth hemispherical, elliptical and 
toriconical heads skirt is no more the -design criterion. 
Fatigue life for vessels v/ith hemispherical and elliptical 
heads become infinite. Fa.tigue life for vessels with tori- 
spherical and toriconical head increase by Sofj and 45 ^ res- 
pectively. It is felt that fatigue life of tori spherical 
head can be further improved by shifting the skirt and modi- 
fying the heating rate of the skirt. 



COHCLUSIOZS 


In this thesis thcrnal fatigue life of vessels v/ith 
foiir t^'pec of heads mmely hemispherical, elliptical, 
torispherical an'd tericonical heads have "been studied# 

Insides of the vessels are s'ahjected to heating, holding 

1 

and cooling and outsides insulated# 

It is observed that mcaciiiiu;.! stress variation occurs 
in the s3cirt region near the head and is the dominant 
factor for the therml fatigue life of the vessels with 
all heads. Of tlie four types of vessels studied, the 
vessel v/ith heaispherica,! heads has the highest thermal 
fatigue life. The thermal fatigue life of the vessels) 
dccj’c.asoc in the or do: of clUiptical, torispherical and 
toriconical heads. 

If stress variation only in iiie heads is considered, 
it is again observed tlmit hemispherical heads have the 
minimu3ii variation. 

It is seen thit thermal fatigue life decreases with 
increase in heating/cooling rate and hold time as expected. 

The decrease is least in the vessel v/ith hemispherical heads. 
The thermal fatigue life decreases in the order of elliptical, 
torispherical and toriconical heads. 

Two methods were studied to reduce iiie variation of 
thermal stresses in tlie skirt : near the head. One was 
changing the position of the skirt and other was to subject 
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the one fourth of the shirt adjacent to the head to the 
same temperature variation ac tlmt of the inside of the 
vessel. She first method does not reduce stress variation 
much. The second method redticos variation considSL’ahly. 
This reduction is dramatic in case of spherical and 
elliptical heads. 
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able 5: Maxiiami Principal Stress, MPa 





























Haxini'um Stress Amplitudes (MPa) in Various 
Pressure Vessels with heating/cooling rate 






































Table 5i Maximton Stress Amplitudes (I'lPa) ia Various 
Vessel Heads with a heating and aooling 
rate of 1o5®C/mino 
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Table 6: Maximm Stress Amplitiide and Fatigue Life 
cf pressure Vessels witli Variois Heads 
for Different hold periods. 
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Table 7; I-Iaximxiia Stress Amplitude (HPa) without and with 
Skirt heating (Heating and Cooling for 4 hourff 
at 1®G/min., hold period nil, temperature 
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Fig. 3 Pressure vessel with ellipsoidal head . 











Fig. 5 Pressure vessel with toriconical head. 
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Thickness, mm 

Fig. 7 Maximum principal stress across A A in upper heads due to 
mechanical loading. 




Fig. 8 Temperature distribution at the end of heating for 
vessel with hemispherical head. 
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Fig. 9 Temperature distribution at the end of cooling for the 
vessel with hemispherical head. 
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Fig. 10 Maximum principal stress across AA In the 
upper heads due to mechanical and thermo! 
loading. (Typical Thermal Cycle) 
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Fig. 11 Maximum Principal stress across AA in the 
upper heads due to mechanical and thermal 
loading (Typical Thermal Cycle) 
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Fig. 12 Maximum principal stress across BB in the 
cylindrical portion due to mechanical and 
thermal loading (Typical Thermal Cycle) 
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Fig. 13 Stress cycle at a point in the heads giving maximum stress 
amplitude. 
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Fig. 16 Maximum stress amplitude vs hold period 
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